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^ ! Abstract 

en 

We consider the dielectric breakdown model in the limit 7/ — > 0^. A differ- 
^^ I ential equation describing the surface growth is derived; this equation is KPZ 

c/3 ' plus a term causing linear instability, and includes a short-distance regular- 



ization similar to a viscosity. The equation exhibits an interesting dynamics 
in terms of poles, permitting us to derive a large family of solutions to the 
equation of motion. In terms of poles, we find all stationary configurations 
of the surface, and analytically calculate their stability. For each value of the 



a 

O I viscosity, only one stable configuration is found, but this configuration is non- 

linearly unstable in the presence of an exponentially small amount of noise. 

■ The present approach may be useful in understanding the dynamics for finite 

^ , r/, in particular the DLA model, in terms of perturbations to the infinitesimal 

'nT I r] problem. 

(N 
O 

I. INTRODUCTION 

G\ 

The problem of Laplacian growth in two dimensions has been much studied since the 
ci I introduction of the model of Diffusion Limited Aggregation (DLA) [jl| . The model is known 

to give rise to complex, branching structures, with a non-trivial fractal dimension. 

rj^ • In radial geometry, the DLA model is defined by starting with a seed particle located in 

C , the two-dimensional plane, and letting additional particles arrive from infinity by a series of 

Q ! random walks. Equivalently, one can define a Laplacian field, 0, subject to V^0 = 0, with 

boundary conditions such that = on the cluster boundary and 0(r) ex log(r) as r — > cx). 

(,j Then, the DLA model is defined by the statement that the probability of adding a particle 

'j_j ■ at any point along the cluster boundary is proportional to the normal derivative of 0, where 

one solves for after every growth step. 

One can also define the problem in strip geometry, which is what we will use throughout 
this paper. Here, the growth is confined to a strip, which is periodic in the horizontal 
direction, with the growth happening vertically. The initial, seed cluster is taken to be a 
flat surface, and random walkers are released from infinity, above the strip. In the strip 
geometry, one finds fractal clusters when the ratio of the width of the strip to the walker 
size diverges. We will refer to this as the IR limit of the problem. 

A useful generalization of the DLA model is the Dielectric Breakdown Model (DBM) , 
which includes an additional free parameter 77. Then, the probability to add a particle at a 
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given point along the cluster boundary is taken to be the normal derivative at that point, 
raised to the power t]. For large 77, growth occurs largely at the tips of the cluster, and the 
branches of the cluster become very long and thin, with a fractal dimension approaching 1 as 
1] -^ 00. For small 77, the cluster becomes less branched, with a dimensionality approaching 
2 in the limit rj —>■ 0, where the dynamics reduces to an Eden model |^. 

The idea behind the present work is to examine the DBM in the limit of infinitesimal, 
positive rj. For any positive rj, in strip geometry, the flat surface is unstable. However, in 
the limit r^ — > 0^, the fluctuations about the fiat surface become small. In the next section, 
we will show how to rescale the fluctuations in the surface to obtain a non-trivial growth 
equation describing the rj ^ 0^ limit. 

What we hope to do is to describe the rj ^ 0^ limit sufficiently well that the problem 
at finite r] can be handled by some kind of perturbation theory; we are still far from this 
goal, but this paper represents a first step in that direction. Note, though, that for any 
non-infinitesimal 77, the growth model gives rise to fractal structures with nontrivial scaling; 
this implies that the perturbation theory will give rise to an IR divergence and require RG 
techniques. We will discuss this more in the conclusion. The nontrivial scaling at finite 77 
should be associated with a crossover length, such that on scales less than this length the 
surface would be self-affine, like KPZ, but on scales greater than this length it would be 
fractal. 

We will find that our problem exhibits an interesting pole dynamics. We derive a large 
family of solutions to the equation of motion using the pole dynamics, and find all stationary 
configurations of the surface, as well the stability of these configurations. We will find that 
the 7/ — >■ 0^ limit of the problems contains many of the essential physical features of the 
finite t] problem: non-linear instability for exponentially small bare noise and competition 
between fingers. Finally, the dynamics of moving poles that we will uncover is very similar 
to the branched growth dynamics of Halsey and coworkers |Q ; the location of the poles that 
we find will correspond to points on the surface of the cluster separating branches in that 
approach. In a future work |^, we will elaborate more on the relationship of pole dynamics 
to branches of the cluster, and present a more complete theory of the statistical properties 
of the surface when the system is destabilized by noise. 

In addition to being a useful limit of the DBM, the 77 ^ 0+ model should be useful for 
certain other physical situations. The instability term that we will derive is to some extent 
reminiscent of a MuUins-Sekerka instability 0, found in solidification problems. While the 
DLA model describes growth processes limited by a diffusive field, the KPZ model describes 
growth of a surface normal to itself, with no diffusive field. One expects that the equations 
we derive will be valid in describing growth of a surface normal to itself, with some weak 
coupling to a diffusive field, such as would be found in solidification problems. 

II. DIFFERENTIAL EQUATION FOR SURFACE GROWTH 

In this section, we derive the differential equation describing the surface growth in the 
?7 — > limit. First we introduce the continuum equation of motion for the DBM, and then 
we take the appropriate limit. We rescale the fluctuations in the height of the surface in 
this limit. We regularize the equation on a short distance, and argue that this regulariza- 



tion corresponds to a finite walker size in the DBM. Finally we discuss some basic scaling 
properties of the resulting equation and discuss how to include noise in the dynamics. 

Although some time will be spent deriving the equation of motion in the rj ^ limit, 
the final result is quite simple. It will turn out to just be KPZ dynamics 0], describing the 
surface at r/ = 0, with an additional linear instability due to the non-vanishing rj. 

Consider the dielectric breakdown model. As discussed above, there are two different 
kinds of geometries commonly used. One possibility is to grow out radially. Another possi- 
bility is to grow upwards on a strip of finite width, using periodic boundary conditions. We 
will use the strip geometry throughout this paper, although in later sections we will occa- 
sionally discuss limiting cases of the strip geometry, in which the width of the strip becomes 
large, and the growth becomes the same as growth upwards from the real line. For the strip 
boundary conditions, we will use a coordinate x, periodic with period 2it, to parametrize 
the cluster. 

We will use complex variables to describe the coordinates in the two-dimensional strip 
in which growth occurs. The real part of the complex variable will represent the horizontal 
position, while the imaginary part represents the vertical position. Growth will occur in the 
positive imaginary direction. 

We will parametrize the boundary of the cluster by using a function F{x), such that F{x) 
yields the cluster boundary, for x = 0...27r. F{x) is taken to be analytic, and one-to-one 
for X with positive imaginary part. In the limit that x -^ +ioo, we require that F{x) -^ x. 
This implies that F{x) can be written as 

oo 

F{x)=x + J2Fij)e'^^ (1) 

j=0 

Throughout, we will use the symbols j, k, ... to refer to modes in Fourier space, while x will 
be used to refer to real space. The function F{j) vanishes for j < 0. 

The surface does not contain overhangs in our limit, although overhangs are essential for 
the fractal clusters at finite rj. The absence of overhangs permits us to describe the surface 
by its height, a function h{x). If we expand h{x) in powers of F we obtain 

h{x) = Im(F(x)) + ... (2) 

For small F, the higher order terms can be neglected. Similarly, for small h, we find that F 
is just the positive wavevector part of h. At finite t], the equations of motion in terms of h 
will become singular due to the presence of overhangs, but the equation of motion in terms 
of F will not have any such problems. 

At given rj, the dynamics of the cluster can be obtained from the Shraiman-Bensimon 
equation 0]: 

d,F{x,t) = z{d.F{x,t))Jdx'\d.,F{x')\-'-'^^^r^ (3) 

Let us introduce a regularization into the equation. We would like to add a term that 
cuts the problem off at a constant length in real space, corresponding to the finite walker 
size in a DLA or DBM problem. However, for small F{j), the constant cutoff in real space 



can be replaced by a constant cutoff in x space. A constant cutoff in x space can be obtained 
by writing 
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dtF{x, t) = t{d.F{x, t)) J dx' \d.>F{x' + iw)\''-'' l,,,_l,, (4) 

for some u. While equation (|]) leads to finite time singularities, equation @ does not. 

Let us write this equation in terms of F{j) and expand this equation in powers of rj and 
Ofc. The linear term is 

d^F{j) = r]kF{k) - r]uk'^F{j) (5) 

Unlike the linear term, the nonlinear terms in equation (^) are nonvanishing in the limit 
r; — > 0. The first nonlinear term is, at ?] = 0, 

dtF{j) = -J2 {Hj - k)F{k)F{j -k)- 2k{j + k)F{k)F{j + k)) (6) 

A;>0 

There is also a constant term in the expansion of equation (^ which simply describes an 
overall upward motion of the surface. We will drop this term by going to the comoving 
frame of the surface. 

For small r], so long as the combination of equations (|],^ leads to well behaved solutions, 
we expect that F{k) will be of order t], simply by scaling. This implies that all higher order 
non-linearities will be unimportant in the limit. Since F{k) is taken to be of order rj, the 
fluctuations about the fiat surface vanish in the limit of vanishing rj, and so the growing 
cluster has fractal dimension 2 in this limit, as explained in the introduction. 

So, if we rescale the field F by r], and also rescale the time coordinate, the equation of 
motion for the system will be 

dtFU) = kF{k) - vk^F{j) - Y: {Hj - k)F{k)F{j - k) - 2k{j + k)F{k)F{j + k)) (7) 

A;>0 

While the full equation (^ leads to dynamics without finite-time singularities, we have 
not shown this for equation (0) obtained from equation (|). However, all of the solutions we 
write down are well-behaved for all time, so we will not worry about this point too much. 

Using equation (0) and equation (|^ we can obtain a differential equation describing the 
growth of the height. The result is 

dth{x, t) = \d^\h{x, t) + udlh{x, t) + {d^h{x, t)f (8) 

This is equivalent to equation (|^). This equation is similar to a KPZ equation 0, except for 
the presence of the term \dx\h{x, t), which leads to an instability. The operator \dx\ is a non- 
local operator in real space, equal to \j\ in Fourier space. The presence of this term change 
the dynamics greatly from the KPZ problem. In the KPZ problem, all states eventually 
relax to the state h{x,t) = const, in the long time, zero noise, limit. The instability implies 
that the state with constant h{x) is not stable to small perturbations. 

Let us discuss the possibility of perturbing equation (H) with noise. We can easily modify 
the equation to 



dth{x,t) = \d.^\h{x,t) + iydlh{x,t) + {d^h{x,t)f + noise (9) 

where some random noise field has been added. We expect that the noise will be short range 
correlated on a length and time scale set by the short-distance regularization. However, the 
magnitude of the noise should also depend on the short distance length scale. The physical 
origin of the noise is shot noise due to fluctuations in the number of walkers arriving at a 
given point; the size of these fluctuations will depend on the size of the walkers, as for a 
large number of walkers the shot noise will get averaged out. To get at the magnitude of 
the noise, let us discuss some scahng properties of equation (H). 

We assumed from the beginning boundary conditions such that x was periodic with 
period 27r. To discuss scaling, we will change these boundary conditions; we will now 
measure x in terms of a short- distance length scale set by the regularization, so that x is 
periodic with period 27r/i/. The natural differential equation to describe the dynamics at 
short distances, of order the walker size, is 

dth{x, t) = \d,\h{x, t) + dlh{x, t) + {d^h{x, t)f + N{x, t) (10) 

where the correlation length and time of the noise field N{x, t) are now of order unity, 
and the mean-square fluctuations in N{x, t) are also of order unity. The magnitude of 
the fluctuations will be set by noise- reduction parameters 0, but will not depend on the 
macroscopic length scale, 2tt/v. Now, let us rescale length and time by a factor of v, while 
rescaling h{x,t) by a factor of l/u. The result is 

dth{x, t) = \d^\h{x, t) + vdlh{x, t) + {d^h{x, t)f + iV(x, t) (11) 



where N is short-range correlated with mean-square fluctuations of order z/^. Equation (|TT]) 
is the same as equation (H), with an appropriate value of noise added. As v gets smaller, 
meaning that the ratio of the IR cutoff to the UV cutoff gets larger, implying that we are 
approaching the true long-distance behavior, we see that the terms i'd'^h{x,t) and N{x,t) 
get smaller, suggesting that these terms could be dropped in the IR limit. However, the 
term ud1h{x,t) is dangerously irrelevant, so it would be too quick to assume that this term 
(or the noise term) could be dropped; in fact, we will find later that the noise becomes very 
important in the IR limit. 

III. POLE DYNAMICS AT NON- VANISHING u 

We will look for solutions to equation (^. The natural way to describe these solutions 
is in terms of the poles of dxF{x,t). First we will find a solution with one pole and then 
describe the dynamics of a solution with several poles. This dynamics is similar to Calogero- 
Sutherland dynamics, but different enough that the dynamics of the system of poles is not 
integrable. At the end of this section, we also consider the system with different boundary 
conditions, so that F{x, t) is defined for all x on the real line, instead of just on a strip of 
width 27r. The pole dynamics takes a slightly simpler form in that case. 

Similar pole decompositions have been found before, in the KPZ equation [|T0| and KdV 
equation [|ll|. The pole decomposition in continuum DLA is also well known |8|JT^. The 



interest in the particular pole decomposition here is that the equation of motion we are 
considering is regularized (unlike that in continuum DLA) and is unstable when considered 
near h{x) = (unlike KPZ). 

First let us give a simple solution to the equation of motion. From the scaling arguments 
at the end of the last section, we can guess that we want F{k) to behave as 1/k for small k. 
The solution we will pick is 

F{k) = -%-''^ (12) 

rC 

for k > 0, with e being negative imaginary. One may verify that, if equation (^) is inserted 
into equation (|^, the result is 

g2Im(e) 

d,e = z-tiy- 2tiy^-—^j^ (13) 

We find that if e = — ilog(^^4ii^) then dte = and we have a stationary solution to the 
equation of motion. In the small z/ limit, e ~ —iu. 
From equation (|l^), we find that 

p—ie+ix 

d^F(x) = -iu- r-— (14) 

In the small e limit this reduces to 

1 



— zu- 



(15) 



te — IX 
So, dxF[x) has a pole at x = ie. Also we find that 

F{x) = z/log(l - e-"+") (16) 

h{x) = -log(l + e^i"^^^) - 2e^'"(^)Cos(2; - Re(e))) (17) 

Therefore, for x near Re(e), the height h{x) has a very deep depression. The function h{x) 
is negative and diverging logarithmically (up to a cutoff set by Ime) as x approaches Re(e). 
For X near Re(e) + it, the function h{x) has a very broad finger. One should expect that 
the fingers would be broad in this problem, as rj is taken to be infinitesimal and we know 
physically that at smaller rj the fingers get wider. A picture of h{x) for this single pole 
solution is shown in figure 1. 

Under the dynamics, the position of the logarithm in h moves, but its magnitude does 
not. This is interesting, since it suggests that e'^^^^^'^, a quantity that arises in a Cole-Hopf 
transformation of Burgers equation, has simple zeroes. 

We can consider a more complicated solution combining A^ different poles. Consider 

^ u 
F{k) = - ^ -e~'^"^ (18) 

a=l ^ 



In this case, the equation of motion yields 

2Im(e„) 1 1 

d^ea = i-iu-2iu( ^TT^^ + Yi ■■ — - + - = )) (19) 

The dynamics of equation (|19]) causes poles to attract in the real direction while repelling 
in the imaginary direction. Consider a problem with only two poles. If the two poles 
are nearby, with the same imaginary coordinate, but different real coordinates, the poles 
approach and collide in a finite time. However, this singularity causes no problem: the 
poles hit, and then move away from each other in the imaginary direction. Until the poles 
collide, both poles have the same imaginary coordinate and different real coordinates. After 
collision, they have the same real coordinate and different imaginary coordinates. If the 
initial imaginary coordinates are even infinitesimally different, the poles never collide. The 
motion of this two pole problem is sketched in figure 2. 

The physical mechanism for the attraction of the poles in the real direction is simple. 
As discussed above, a given pole causes a large depression in the surface. The surface slopes 
down towards this depression on both sides. This gradient in h then attract the other poles. 

It is interesting to note that the attraction of poles in the real direction provides a 
physical interpretation for the fact that each pole contributes an amount of order v to h{x) 
(see equation (0). Each pole can also be thought of as having a "self-attraction". The term 
d^hix) tends to smear out poles in the real direction; only if a pole has weight at least v 
can it attract itself enough to avoid being spread out. The term in the nonlinear equation 
of motion which causes the interaction between poles is exactly the same as the term that 
cancels d'^h{x) in the single pole dynamics. 

Without the term —i + iu in equation (|I^, the equation could be reduced to a com- 
plexification of Calogero- Sutherland dynamics [|l3,|l^, and then the pole dynamics would be 
integrable. It seems to us that the term —i + iv destroys integrability in equation ([TP|). 

The dynamics of equation (|19D assumes a simpler form if we assume that the poles are 
all very near to each other with e and v small, or, equivalently, if we change boundary 
conditions so that x is now defined on the real line and is not periodic. Then we find 

dtea = t - MttTT-^ + E(^— + -^)) (20) 



IV. STATIONARY STATES OF POLES 

In this section, we consider the qualitative features of the pole dynamics found in the 
previous section. First, we consider the case with F{x,t) defined on the real line, and find 
the long time limit of the pole dynamics, for any finite number of poles. It is found that, for 
given number of poles, the long time limit of F{x, t) is independent of the initial positions of 
the poles, up to a trivial translation in x. Then, we consider the case with periodic boundary 
conditions, and find that, for given initial number of poles, the long time limit of F{x,t) 
tends to one of a discrete number of different states, again up to translations in x. It is 



argued that all except for one of these states is unstable to small perturbations in the initial 
configuration of poles. 

As discussed in the previous section, the poles attract each other in the real direction, 
and repel in the imaginary direction. If we have a finite number of poles, with F{x, t) 
defined on the real line, it is clear that, in the long time limit, all poles will have the same 
real coordinate (one can easily show this from the equation of motion by taking a time 
derivative of I]a(R'6(ea))^). Although the poles repel each other in the imaginary direction, 
they are attracted to the real axis. The result is that the poles form some stable configuration 
with all poles having the same real coordinate, and with different imaginary coordinates. 
There is only one stationary configuration for a given number of poles. This configuration 
is stable in that any small perturbation of the pole positions will decay in time, except for 
the trivial zero mode associated with translation. If this configuration has A^ poles, then 
h{x) behaves as i>Nlog{x) for large x, so configurations with different numbers of poles have 
different boundary conditions at large x. 

Consider stationary states of poles with periodic boundary conditions. One possibility is 
for all poles to have the same real coordinate, up to multiples of 27r. There are, however, other 
possibilities. One very simple possibility is to take 2 poles, with Re(ei) = and Re(e2) = vr. 
In this case, the poles can find a stationary configuration, but this configuration will not 
be stable: any small perturbation in the real part of one of the e^ will cause the two poles 
to move towards each other, until a new, stable, configuration is found with both poles 
having the same real part of e. This process is exactly the process of competition between 
two fingers, with one finger eventually winning out, commonly seen in DLA and Hele-Shaw 
problems. We show some stages of this process in figures 3 and 4. There are of course more 
complicated possibilities for stationary states, but the only stable states are those with all 
poles having the same real coordinate. 

Let us show for strip boundary conditions that there is a maximum number of poles one 
can use to form these stable configurations. Consider what equation (^^ implies for such a 
configuration. Clearly, Re{dtea) = for all a. We also need lm.{dtea) = for all a. Let us 
order the poles so that if a < 6 then Im(ea) > Im(eb). Consider the equation of motion for 
eAT. We find that 

„2Im(eiv) 1 1 

d^eN = i-iiy-i2u( ^^ + y( ., — — + - ^ )) (21) 

The first term on the right-hand side is positive imaginary, while all the other terms are 
negative imaginary. So long as Im(eAr) < Im(ea) for all a ^ N, then the imaginary part of 
the right-hand side is monotonically increasing as Im(eAr) becomes more and more negative. 
In the limit cn — ^ —ioo, the right-hand side becomes 

i{l - u - 2u{N - 1)) (22) 

If this quantity is negative imaginary, then it is not possible, for any cat, to have dte^ = 0. 
Then, 9tlm(e7v) < for all time, and we find that e^ moves off to —ioo. In the limit that 
ejsf moves off to infinity, this pole no longer contributes to the sum in equation ([T8|) and can 
be ignored. 



Putting all these arguments together, we find that, for any v in the strip geometry, we 
can find states with A^ poles, stable against any small perturbation in the pole position. For 
each A^, there is only one such state, and it has all poles with the same real coordinate. 
Such a stable solution with A^ poles can only be constructed if — 1 + z/ + 2v{N — 1) < 0. 
Otherwise, some number of poles move off to infinity until the solution reduces to one with 
-l + z/ + 2z/(iV-l) <0. 

V. STABILITY ANALYSIS 

We will look at the linear stability of the stationary states found in the last section, for the 
case of periodic boundary conditions. Here, we will consider arbitrary small perturbations 
of /i, including perturbations that cannot be written in terms of a small change in the pole 
coordinates. We are able to analytically calculate the eigenvalues of the linear problem, and 
find that, for each value of viscosity, only one of the stationary states is linearly stable. After 
obtaining the result, it is argued that the qualitative features of the result could be obtained 
from the analysis of the stationary states in the previous section: the linearly stable state is 
the only state stable towards adding a pole. 

Before doing any mathematical analysis of our problem, let us recall the stability analy- 



sis of the Saffman- Taylor system [Q. Saffman and Taylor found a one-parameter family of 
stationary solutions to the zero surface-tension problem. Surface tension acts as a singular 
perturbation to this problem. Eventually it was shown |T5| that a small amount of surface 
tension stabilizes the "1/2" finger against small perturbations. However, in the DLA prob- 
lem, although the system is regularized due to non-vanishing walker size, a fractal structure 
emerges instead of stable fingers. In our system, although the system is stable, in the IR 
limit an exponentially small amount of noise will destabilize it, as will be discussed more 
in section VII. In that case, it may turn out that we obtain a complicated structure in our 
problem, for a small amount of noise. 

The physical reason for stability is most easily understood in a WKB approximation. 
Consider a locahzed, short- wavelength, perturbation f{x) to a stationary state, Fq{x). Let 
the stationary state be such that all poles have imaginary coordinate 0. Then, the surface 
has a deep depression near a; = 0, and a very broad finger located around x = ■n. The 
perturbation will move in position, due to the term in the equation of motion dxf{x)dxFQ{x), 
until the perturbation moves to a; = 0, where it disappears into the deep depression on the 
surface. A similar physical mechanism for stability is known in the Hele-Shaw problem and 
other problems [1^. The closer the perturbation is to the tip of the finger initially, the 
longer it will take it to drift along the side of the finger and disappear. Fortunately, the 
non-vanishing viscosity prevents us from localizing a perturbation exactly at the tip as the 
perturbation must have a width of order u or greater if it is to be unstable, so eventually 
all short-wavelength perturbations will be destroyed. 

Now, let us proceed to a more careful analysis. Let Fq{x) be a stationary state, and let 
us consider a small perturbation f{x). We will consider the evolution of the system in time, 
to linear order in /. One possibility to do this is, of course, to take the equation of motion 
(0), writing F{x) = Fq{x) + f{x), and directly derive the equation of motion to linear order 
in /, obtaining dtf{x,t) = Lf, where L is some linear operator whose eigenvalues describe 



the stability of the state Fq. However, this procedure would be very awkward, and we will 
use a more clever technique. 

Let Fq{x) be a state formed using A^ poles, at positions ei, 62, ..., eAr. We will write a 
state F{j), near to Fq{j), as 

f(i) ^ • . 

F{j) = FoU) + ^-^-iY: z//„e-»^ (23) 

3 a=l 

where f{j) is a function of j and the various fa are numbers representing small deviations of 
the pole coordinates from their original positions. By doing this, we have introduced some 
redundancy in describing the possible perturbations to Fq. 

If /(j) = 0, then we can derive a linear equation of motion for the fa-, simply using the 
pole dynamics equation. We must place the N poles at positions ei + /i, €2 + /2, •.., Cat + /at 
and linearize equation ([T9|) about /a = 0. We will obtain some equation of the form 

dtfa = La,tfb + Ma,J, (24) 

for some linear operators L, M, with an implied sum over b. The result of this is that there 
are 2N different modes (there are N complex coordinates fa). One of these modes is a zero 
mode; this is a mode with all fa equal to the same real constant. From the analysis of last 
section, if all e^ have the same imaginary coordinate, then the other 2N — 1 modes are all 
stable (have negative eigenvalue). 

Now, consider the case with non-zero f{j). For simplicity, first consider a case with 
initial conditions such that /(I) is non-zero, but all other f{j) are zero. Then, the linearized 
equation of motion yields 

N 

d,{F{j) - FoU)) = (1 - iyiy,,if{l) + 2z/^((l - 5,,i)/(l)e-^(^-^) + /(l)e-'(^+^)) (25) 

The reason for the factor of {l—6j^i) in the sum is that in equation ([T^), each pole contributes 
only to F{k) for k > 0. It is convenient now to rewrite equation (pSf) as 

N 

dt{F{j) - FoU)) = (1 - ^1' - 2i^N)5,,J{l) + 2z.^(/(l)e-'(^-i) + /(l)e-^^»(^+^)) (26) 

which, combined with equation (P^, is equivalent to the equations 

dJil) = (1 - ul^ - 2uN)f{l) (27) 



dtfa = -2t{f{l)e''^ + f{l)e-''^) + La,bfb + M,, J, (28) 

If we extend this procedure to the case of general /(j), we find 

N 

dtfU) = (j - i^f - 2juN)fij) + 2^^.^ ^(e^(^-^)- - e-'^'~^^'-) f (k) (29) 

k>j a=l 
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dtfa = -2i J2ifij)e'^'^ + fU)e-'^'^) + L,,th + M^^.f, (30) 

j 

The above two equations fully define the linear evolution of the system. Note that 
the matrix describing the linear evolution of the system is triangular, we find that dtf{j) 
depends only on f{k) for k > j. This makes it possible to directly read off the eigenvalues of 
the matrix. There are 2N eigenvalues which are just the eigenvalues from the evolution of 
equation (|2^) . Then there are eigenvalues which are j{l — 2vN) — vj'^ with j = 1, 2, 3, ... Each 



of these eigenvalues must in fact be counted twice, since there is one such eigenvalue with 
/(i) pure real, and one such with f{j) pure imaginary. So long as we consider the stationary 
state found in the last section, with N poles, all with the same imaginary coordinate and 
with N the largest integer less than ^^, then all eigenvalues of the linear evolution are 
negative, with one (or three) exceptions. The first exception is the zero mode corresponding 
to changing the imaginary coordinate of all the poles by the same amount. The other two 
possible exceptions occur if ^^ is an integer, in which case 1 — 2i'N — ul'^ = 0, and so we 
have two more zero modes. However, in any case, there are no positive eigenvalues. 

The above analysis might be slightly confusing, as it seems we have introduced more 
eigenvectors than we started with by adding the fa coordinates. If we look at the linearized 
equation of motion in its original form, without introducing the additional coordinates fa, 
we notice that the diagonal term in the linear equation of motion is j — uj"^. For very 
large j, this term must dominate all other terms, and so, if we order the eigenvalues of the 
linearized equation of motion, and look at the 2j-th eigenvalue, this eigenvalue must be 
close to j — z/j^. This gives us one way to count the number of eigenvalues of the problem: 
since the dimension of the space of f{j) is infinite, we cannot simply count the number of 
eigenvalues directly, but we can count the number of eigenvalues less than a given number. 
Now, notice that 

j{l~2iyN)-uf = k-uk'^ -N + N^u (31) 

with 

k=j + N (32) 

The quantity —N + N'^u will be of order A^, but for very large k it will be small compared to 
k — uk"^. So, while we have introduced 2N eigenvalues by adding the fa coordinates, equation 



(B3) makes it clear that we have made up for this by losing 2N eigenvalues elsewhere: the 



full set of eigenvalues includes the 2N eigenvalues from equation (|2^) as well as eigenvalues 



j - uf - N + N^u with j = N+l,N + 2,N + 3,.... Then, if we look at the 2j-th eigenvalue, 
it will be close to j — z/j^. It is possible to make this argument more precise to show that 
we have found all eigenvalues of the linear evolution. 

The linearly stable configuration is also stable to adding a pole. If a pole is added, one 
pole will move off to —ioo, and the solution will be unchanged. This provides some insight 
into the reason for the linear stability of the solution. One way to create an infinitesimal 
perturbation to h using only poles is to add an additional pole at a very large negative 
imaginary coordinate. If N is large enough, the added pole will move to infinity; otherwise, 
the added pole will approach the real axis and F will not be stable to this perturbation. 
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VI. FLUID DYNAMICS 



Consider the limit of small z/. From equation (|TB|), each pole contributes an amount of 
order v to F{x). However, from the analysis of the previous two section, we expect that the 
only stable state is a state with A^ poles, where, for small v^ we find ^ ~ ^- So, in the 
limit of small z/, we will find that Fix) tends to a finite, non- vanishing function. We expect 
that, in the limit of vanishing z/, the problem can be described in terms of a "fluid" of poles. 
Let us analytically find the state of this fluid which corresponds to the small v limit of the 
stable state found in the last section. 

For finite v and A^, taking all poles at the same real coordinate, the equation of motion 
(p!9|) yields A^ different equations which constrain the imaginary coordinates of the poles. 
These equations are non-linear, and not easy to solve. However, for small v^ the number of 
poles becomes large, while the separation of poles, and the contribution of each pole to the 
right-hand side of equation (|T9|) , becomes small. In that case, we can rewrite (|l^) in terms 
of a density of poles. 

First, let us write the equation for the pole density in general, for any state, not just a 
stationary state. Let the pole density be a function p(e), defined such that 



F{k)= I dhp{e)e-'''' (33) 

We find that the equation of motion ([T9| ) becomes 

d.pie) = -Re(d,{p{e)J{e))) (34) 



where 

Ji^) = ^ - 2^ / d^^' P(e') ( Y^i^ + ;i^^) (35) 

This can be found either as a limit of equation (p!9D, or by inserting equation (^) into the 
zero viscosity limit of equation (^. Here J(e) represents the velocity at which the fluid 
moves, so the equation of motion for the fluid is just a continuity equation. 

In the zero viscosity limit, it is not obvious that h{x) will always be well-defined. For 
finite viscosity, poles cannot hit the axis, but for vanishing viscosity, if the fluid of poles 
moves so that it intersects the real axis in a sufficiently singular manner, equation (|33D will 
become ill-defined. 

Let us look for the stable configuration of poles. Let all poles have real coordinate equal 
to 0, so that p(e) = 5(Re(e))/(— Im(e)). A stationary configuration of the fiuid is one such 
that J(e) vanishes for any e such that p(e) is nonvanishing. For the stable configuration, 
which is the one with the maximum number of poles, p(e) is non-vanishing for all e on the 
negative imaginary axis; that is, there is a density of poles extending off to infinity. 

So, to find the stable configuration, we must find a function / such that J vanishes on 
the negative imaginary axis. That is, we must solve the integral equation 

-l + 2/d,/(y)(^— j^ + -j^;^)=0 (36) 
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for all real x > 0. Equation (^) can be rewritten as 

7 1 _ 1 







Finally, introducing a new coordinate, z = Cosh(x), and z' = Cosh(y), we find 



oo 



V7^^'^''l^^7^ '^«' 



for all 2; > 1. This is an integral equation which can be solved using known techniques for 
singular integral equations [ p!7|] . 



We find that 

^, , 1 , ,l + Cosh(x), , ^. 

This function has vanishes exponentially for large x. This implies that the density of poles 
vanishes exponentially for large, negative imaginary coordinates. 

For this given /, one can examine the function h{x) that results. As x — >■ one finds 
that dxh{x) diverges logarithmically. However, the height h{x) does not diverge as a; — >■ 0. 
It is interesting that the fingering solution in the zero viscosity limit is so well behaved for 
?7 ^ 0^; for r] = 1 it is known that the stable finger only occupies half the width of the 
channel, going off to infinity at two distinct points. It is possible that, for any non-zero rj, 
the solution we have found will narrow so that it does not occupy the whole channel. 



VII. NOISE IN THE IR LIMIT AND PERTURBATION THEORY 

From the naive scaling analysis of equation ( |lO|) , we concluded in section II that the noise 
vanishes in the IR limit. However, as we noted, the viscosity also vanishes in this limit, and 
since the viscosity is a dangerously irrelevant operator, this can change everything. In this 
section we will argue that, in fact, in the IR limit, the system becomes unstable to an 
exponentially small amount of noise. Then, we will briefiy discuss the perturbation theory 
needed to deal with finite 77. However, the lack of stability of the 77 — >■ O"*" system in the IR 
limit implies that we should not apply the perturbation theory in 77 until we have first found 
some way to deal with noise at infinitesimal 77. 



Consider the equation of motion ([T0|), so that we are using the short- distance length scale. 
On this length scale, the stationary configuration of h{x), found above, appears to be very 
smooth. So, for the interaction of the stationary configuration with the noise, we can use a 
WKB approximation, as discussed above. In the WKB approximation, the perturbations to 
the surface move along the surface until they hit the depression caused by the poles, where 
the perturbations are destroyed. 

In the rescaled lengths, it takes a time of order v~^ to hit the depression. So, any small 
perturbation will have a long time to grow before it is destroyed. In rescaled lengths and 
times, the most unstable wavevector, when considering fiuctuations about a fiat surface, is 
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j = 1/2, which grows at a rate equal to 1/4. So, a short distance perturbation to the surface 
will grow at a rate of order 1/4, for a time of order 1/z/, before being destroyed. So, the 
perturbation will be amplified by a factor such as e^^. So, an exponentially small amount 
of noise will lead to non-linear instabilities of the surface. 

As discussed in section V, it is impossible to localize a perturbation exactly at the tip 
of the finger. However, the greatest growth will occur for perturbations located at the tip 
of the finger, a distance of n away from the depression, as these perturbations will last the 
longest before being swept into the depression. It is expected that perturbations near the 
tip lead to a tip splitting instability. 

The exponential non-linear instability can also be seen from the linear stability analysis 
above; many of the eigenvectors are almost degenerate. In forming the upper triangular 
matrix used in the linear stability analysis, many of the elements above the diagonal are 
exponentially large. 

Let us briefly discuss the theory at finite rj. One must extend equation (|]) to include 
higher terms in 7] and higher non-linearities in F. This leads to a perturbative correction 
of equation (|^ in powers of rj. The perturbation theory should include three effects. First, 
the stationary configurations found in the last section should change even in the absence 
of noise. Secondly, when we drive the theory by a small amount of noise, we have argued 
above that the configuration becomes unstable in the IR limit. For non-vanishing t] the 
fiuctuations in h{x) will interact differently than at vanishing rj. However, until we can fully 
understand the limit of vanishing rj, there is no reason to consider the perturbation theory 
in detail. Thirdly, we argued that in the limit of vanishing t] is was acceptable to replace 
the cutoff of constant length in real space by a cutoff of constant length in x space, where 
X is used to parametrize the surface. This is not valid for finite r], and must be corrected. 

VIII. CONCLUSION 

In conclusion, we have considered the problem of the dielectric breakdown model in the 
limit of vanishing rj. A dynamics in terms of poles was found, and was used to find the 
stable configuration of the system for different values of u. In the limit of vanishing u, this 
provides us with a fingering solution similar to the fingering solution of the Saffman- Taylor 
model. As expected, the finger is much broader than the Saffman- Taylor finger. 

It was argued that in the IR limit an exponentially small amount of noise destabilizes 
the system. This will be discussed will be discussed more in a future work, where we will 
also elaborate more on the branching properties of the surface. 

Once the t] ^ 0^ problem is solved in the presence of noise, an even more interesting 
and difficult problem is to obtain perturbative results for the DBM in powers of r]. Since 
the problem for nonvanishing t] gives rise to fractal structures, with nontrivial scaling, per- 
turbative results for the DBM should involve further IR divergences, which would have to 
be handled by an RG procedure and would give rise to the non-trivial dimensions. We hope 
that the work in this paper represents a first step towards solving that full problem. 

Note Added: After this paper was first posted to cond-mat, Vincent Hakim pointed out 
to me some closely related results in the literature on the problem of fiame fronts. Equation 
(§) is well known in that problem, where it is referred to as Sivashinksy's equation |jl8[. The 
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pole dynamics was previously identified also [|T9|. It appears, however, that the analytic 



calculation of the linear stability that appears in this paper is new, as well as the connection 
between the exponentially large off-diagonal elements of the stability matrix and the WKB 
stability analysis; some numerical results on the stability of the system are known |2^, and 



agree with the exact analytic calculation here. Further, the connection between the DBM 
and Sivashinsky's equation also appears to be new. 
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FIGURES 




FIG. 1. A plot of h{x) for the trial solution at ly = .2 
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FIG. 2. Sketch of paths of two poles. Horizontal line represents real axis; poles lie in lower 
half-plane. Figure is approximate; actual paths will be rounded out. 
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FIG. 3. Competition between two fingers, short times. Poles have real coordinates .1 and tt — .1. 




FIG. 4. Competition between two fingers, long time limit where one finger is completely sup- 
pressed. 



